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Four fermion processes at e + e~ colliders in a range of energy from 100 GeV to 1 TeV are discussed and 
examples of results obtained with an event generator presented. We also investigate the effects of the inclusion 
of initial state QED corrections and QCD contributions. 



1. Introduction 

Heavy vector boson production will be investi- 
gated at e + e~ colliders in a wide range of ener- 
gies. At LEP II, the mass and the width of the 
charged boson will be measured [1] and searches 
for its anomalous couplings performed [2]. The 
relevant process is, in this case, the production of 



two W's above threshold 

e+ e~ W + W~ . (1) 

At somewhat higher energies other processes like 

e+e" Z Z , (2) 

e+e- W ev e , (3) 

e+e" Z e+ e~ , (4) 

e+e~ — > Z v e v e , (5) 

become important [3]. 



Since the produced heavy bosons immediately 
decay, for all above processes the detected exper- 
imental signal is a four fermion final state, which 
has important consequences. First of all, the 
number of possible final states grows. Secondly, 
the number of Feynman diagrams contributing to 
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a given final state can be very large. Finally, the 
same four fermion final state can originate from 
several signals. For example, all 5 signals (l)-(5) 
can end up with a e + e~v e v e final state. 

Of course, at some energies, a particular set of 
diagrams dominates, and all the others can be 
considered as a background. At LEP II, three di- 
agrams with 2 decaying W's dominate, while, at 
next linear colliders in the TeV regime, diagrams 
for single W or Z production, are responsible for 
the leading effects and the former diagrams give 
small contributions. For an exhaustive treatment 
of all questions concerning signal and background 
see ref. [4]. Here, we only stress that it is not al- 
ways possible to use kinematical cuts to remove 
the background effects. For instance, if two neu- 
trinos are present in the final state, no invari- 
ant mass cuts can be performed and only quan- 
tities summed over all possible neutrino flavours 
are suitable for comparisons with experiment. 

Thus, an explicit calculation of all four fermion 
processes is unavoidable, in order to check the 
effectiveness of cuts to reduce the background and 
to extract physical information whenever these 
cuts are not applicable. 

This rather ambitious program has been car- 
ried out in five papers [4]- [8]. The main result is a 
very fast Monte Carlo program, EXCALIBUR [5,6], 
that, in the limit of vanishing fermion masses, 
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takes into account all possible Feynman diagrams 
contributing to any given four fermion final state. 

Note that neglecting fermion masses implies 
the absence of diagrams where a Higgs boson 
couples to the fermions: but, although we can- 
not compute the Higgs signal in our approach, 
we can at least reliably estimate the background. 
For one diagram in the background one then has 
to introduce a special treatment, i.e. the dia- 
gram where two W's fuse to bb with top quark 
exchange, which has to be taken massive. 

The QCD tree level effects in 4 jets final states 
[7] and the initial state QED radiation [8] are 
also implemented, so we feel confident that our 
program is complete and accurate enough to de- 
scribe, in the framework of the Standard Model, 
the four fermion physics relevant at e + e~ collid- 
ers. 

The outline of this paper is as follows. In sec- 
tion 2 we sketch the method used in computing 
the matrix element and in building the 

Monte Carlo. In section 3, 4 and 5 results 
obtained with the described generator are pre- 
sented. Finally, we summarize our conclusions. 

2. The method 

Since one wants to deal with all possible final 
states, first of all an efficient way to compute the 
matrix elements is needed. 

There are Abelian and non Abelian graphs with 
two distinct topologies [4]. All Feynman diagrams 
can be obtained from those two topologies by per- 
muting the fermion momenta. The number of 
diagrams can be very large (up to 144), so the 
evaluation has to be performed at the level of he- 
licity amplitudes. Using spinorial techniques, one 
single basic function can be found that describes 
the spinorial structure of both Abelian and non 
Abelian graphs. Thus, for massless fermions, each 
helicity amplitude consists of a sum of very sys- 
tematic and relatively compact expressions. This 
is important for the computational speed. 

As for the Monte Carlo integration over the 
final fermion momenta, we have used a multi- 
channel approach. The basic idea is to generate 
the integration variables according to distribu- 
tions that approximately reproduce the peaking 



behaviour of the matrix element squared. This 
reduces the variance of the integrand, and there- 
fore the Monte Carlo error. 

Since the peaking structure of the matrix ele- 
ment squared can be very rich, one set of integra- 
tion variables can only describe a limited number 
of peaks, so that, in general, one has to gener- 
ate each peak with different mappings of random 
numbers, that is different channels. In EXCALIBUR 
those channels are automatically chosen by the 
program, looking at the Feynman diagrams con- 
tributing to the selected process. Also the opti- 
mization of the predetermined probabilities used 
to choose the various channels is automatic, lead- 
ing to an effective increase in the program speed 
by almost an order of magnitude [6] . 

QED corrections are implemented using the 
structure-function formalism [9] . Each of the in- 
coming fermions is assumed to have its energy de- 
graded by the emission of photons parallel to the 
beam. For the energy distribution of the fermion 
after radiation we take a structure function $ that 
incorporates the leading log O (a) and O (a 2 ) ini- 
tial state radiation with exponentiation of the re- 
maining soft-photon effects [8]. Our model for the 
total radiative cross section is then 

l l 

<r(s) = j j dx\ d,X2 <&{x\) $(£2) (To{xiX2S) (6) 


where <7o is the nonradiative cross section and 
x\, X2 represent the energy content of the incom- 
ing fermions after radiative emission. This pro- 
vides an adequate description of the leading QED 
effects. 

When four quarks are present in the final state, 
one has to add the concomitant QCD production 
channels, and also the production of a quark pair 
and two gluons, since both types of final states 
will appear as jets. In our Monte Carlo the for- 
mer contribution is easily implemented without 
additional CPU time by adding gluons wherever 
photons connect quark lines [7] (of course the cor- 
rect QCD coupling and colour structure should be 
taken into account). Finally, the latter process 
can be efficiently computed using the recursion 
relations of ref. [10]. 
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Table 1 

Cross section in pb for e + e~ — ► e + e~ v ellyT v e llyT . Cuts applied: m(e+e~) > 10 GeV, Icos^gil < 0.9, 
E e ± > 20 GeV. 



V~s (GeV) 


150 


175 


200 


WW diagrams only 


.3600 ± .0011 10" 2 


.1181 ± .0002 


.1304 ± .0003 


WW + background 


.4625 ± .0011 10" 2 


.1246 ± .0002 


.1591 ± .0003 


Relative difference (%) 


22 


5 


18 


Table 2 








Cross section in pb for e + e~ 


—> e~v e + 2 jets. 


Cuts applied: E e - q gi > 20 GeV, 


1 COs6' e - ) q> q' \ < 0.9, 


COS (q, q')\ < 0.9, m(qq') > 10 GeV. 






V~s (GeV) 


150 


175 


200 


WW diagrams only 


.1921 ± .0060 10" 1 


.6684 ± .0014 


.7202 ± .0016 


WW + background 


.1513 ± .0042 10" 1 


.6690 ± .0014 


.7490 ± .0018 


Relative difference (%) 


-27 


0.1 
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3. Results without initial state radiation 

In this section we shall present examples of 
numbers obtained with EXCALIBUR. We start 
showing results without QED radiation and QCD 
effects [4]. The actual values for the input param- 
eters are a = 1/129, sin 2 9 W = 0.23, M w = 80.5, 
IV = 2.3, M z = 91.19 and T z = 2.5 (all GeV). 

In the first row of table 1 only the three WW 
signal diagrams are considered, while in the sec- 
ond entry all diagrams leading to a final state 
with e + e~ and two neutrinos are taken into ac- 
count. All neutrino flavours are summed over, 
because, as already discussed in the introduction, 
it is not possible to distinguish them experimen- 
tally. As a consequence, the background for this 
process consists of both an interfering and a non- 
interfering part. The last entry shows the whole 
effect. Usually, one tries to get rid of the back- 
ground contributions by imposing invariant mass 
cuts on the particles that can come from decaying 
W's. Note that this is not possible here because 
two neutrinos are produced. Therefore, any study 
on, say, lepton universality made at LEP II on the 
basis of purely leptonic events like those shown 
in table 1, should take care of the background. A 
last remark is in order. While a large background 
below threshold (22%) is somewhat understand- 
able, the large contribution at 200 GeV (18%) 
may appear surprising. This is due to two com- 
bined effects. On one hand diagrams for single W 



production become important at higher energies. 
On the other hand, at 200 GeV, background Z's 
are produced above threshold. 

In table 2 the semileptonic process e + e~ — ► 
e~v e -\- 2 jets is considered (taking the CKM ma- 
trix unity it is the sum of cross sections with ud 
and cs in the final state). Here, only interfering 
background is possible. Then, since no resonant 
Z's are produced, at 200 GeV, the effect of the 
background is not as large as in the previous case. 
For the process e + e~ — ► pTv^ + 2 jets, where 
also diagrams with single W production are ab- 
sent, the contribution of the background is even 
smaller (0.3%), again at 200 GeV. With semilep- 
tonic final states invariant mass cuts are possible. 
A cut 70 < m(e~p e ), m(qq') < 90 GeV reduces 
the background from 4 % to 1% for the process 
in table 2 at 200 GeV, and only 10% of the events 
are lost. 

Finally, in table 3, the influence of the inter- 
fering electroweak background on a specific four 
quark final state is investigated. 

This concludes our discussion on signal and 
background at LEP II. At higher energies sin- 
gle boson production processes become important 
and the problem is to find a way to disentan- 
gle, for any given final state, the signals (l)-(5). 
In ref. [4] we showed that this can be achieved 
by cutting the invariant masses of the events 
around the vector boson masses. For example, 
at 1 TeV and with the cuts E e - u j > 20 GeV , 
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Table 3 

Cross section in pb for e+e" — > uddu. E {aU partides) > 20 GeV, | cos6( aU pa rtides)\ < 0.9. Moreover 
m(ij) > 10 Gel/ and | cos < 0.9 between all possible final state pairs. 



V~s (GeV) 


150 




175 


200 


WW diagrams only .2141 ± 
WW + background .2317 ± 


.0019 10" 1 
.0010 10" 1 




.7699 ± .0018 
.7745 ± .0018 


.8726 ± .0023 
.8987 ± .0023 


Relative difference (%) 


8 




0.6 


3 


Table 4 

Results on radiatively corrected cross sections 
strategies, for the process e + e~ — ► e~v e ud. The 


and average energy losses, under various calculational 
cross sections are given in picobarns, e and e in GeV. 






= 176 GeV 






strategy <7o 


a 




£ 


£ 


WW,f,l .60111 ± .00032 

WW,s,l,a 

WW,s,l,b 

WW,s,2 

WW,cuts .44651 ± .00092 
all,cuts .45011 ± .00097 


.50490 ± .00032 
.50484 ± .00033 
.50175 ± 00098 
.50258 ± .00097 
.37737 ± .00091 
.37926 ± .00095 


1.162 ± 0.002 
1.172 ± 0.002 
1.167 ± 0.006 
1.178 ± 0.006 
1.192 ± 0.007 
1.149 ± 0.007 


1.175 ± 0.002 
1.170 ± 0.006 
1.181 ± 0.006 
1.195 ± 0.007 
1.152 ± 0.007 



Icos^g- u d\ < 0.9, m(ud) > 30 GeV, the cross 
section for the process e + e~ — ► e~v e ud includ- 
ing all diagrams is 0.0452 pb, while the result 
with only WW diagrams is 0.0143 pb. That is 
a difference of a factor 3. An additional cut 
Mw — 2Tvy < m(e - i7 e ), m(ud) < Mw + re- 
duces both numbers to 0.0079 pb. However, even 
if such a procedure removes all the other signals, 
it should be noted that almost half of the WW 
signal has been cut. Therefore, in order to avoid 
this huge loss in statistics, any study about tiny 
effects (for instance searches for anomalous cou- 
plings) has to be performed directly at the level 
of a complete e + e~ — ► 4 fermions generator. 

4. Results with initial state radiation 

In order to present some salient results we have 
chosen the specific process 

e + e~ — ► e~v e ud , (7) 

which contains the WW pair production signal, 
as well as a non-negligible background [8]. 

Table 4 contains numbers from our Monte 
Carlo obtained under different strategies. The 
quoted values are the nonradiative cross section 
Co, the total cross section with initial state ra- 



diation it, the average 'energy loss' e defined as 
■^(1 — X1X2) (see ref. [11]) and the real energy 

loss ^(2 — x\ — X2), denoted by e. The strategies 
used are 

• WW,f,l: the WW signal diagrams only, 
with the flux- function approach to O (a) as 
given in ref. [11]; no cuts; 

• WW,s,l,a: the WW signal diagrams only, 
with structure functions that are simply the 
flux function of ref. [11] to O (a), in which 
2a is replaced by a; no cuts; 

• WW,s,l,b: the WW signal diagrams only, 
with O (a) structure functions of ref. [8]; no 
cuts; 

• WW,s,2: the WW signal diagrams only, 
with O (a 2 ) structure functions as given in 
ref. [8]; no cuts; 

• WW, cuts: like the previous case, except 
that now we also impose the following cuts: 
E e - > u> s > 20 GeV, \co S e - t u> 3 \ < 0.9, 
I cos (ud)\ < 0.9, m(ud) > 10 GeV; 

• all, cuts: like the previous case, except that 
now also all the background Feynman dia- 
grams are taken into account. 
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Table 5 



Cross section in pb for e + e — ► 


4 jets. 


Cuts as in table 3. 




V~s (GeV) 




175 


190 


WW signal 




3.0674 ± 0.0074 


3.5136 ± 0.0090 


WW signal + ISR 




2.5622 ± 0.0071 


3.1416 ± 0.0089 


All EW diagrams + ISR 




2.5922 ± 0.0075 


3.3553 ± 0.0097 


All EW diagrams + interfering 


QCD -+ 


- ISR 2.6202 ± 0.0079 


3.3789 ± 0.0100 


All EW diagrams + all QCD + 


ISR 


3.1155 ± 0.0123 


3.8688 ± 0.0146 



For this analysis we have adopted the values given 
in ref. [11], namely 

a = (127.29) _1 , sin 2 9 W = 0.2325 , 
a str.f.,flux = (137.036)" 1 , 
M z = 91.173 GeV , T z = 2.4971 GeV , 
M w = 80.220 GeV , IV = 2.033 GeV . 

Referring to eq. 6, a is contained in <7o, while 
a str f flux i s used in <&(x). The structure func- 
tion method differs slightly from the flux function 
method as comparisons between the first and sec- 
ond row of table 4 show. The inclusion of cuts or 
more diagrams affect both cross sections and en- 
ergy losses. 

Since the proposed direct reconstruction 
method for the W mass suffers from a shift in Mw 
due to the radiated energy [1], a precise knowl- 
edge of e and e is warranted. Our results show 
that the precise treatment of ISR, the choice of 
cuts and the neglect of diagrams all affect the en- 
ergy losses. The first effect was also found in ref. 
[11], where also the influence of the Coulomb sin- 
gularity is discussed. The two other effects show 
that a Monte Carlo treatment allowing for cuts 
and being able to include all diagrams is indis- 
pensable. 

5. Results with initial state radiation and 
QCD contributions 

With four-quark final states also the QCD con- 
tributions must be taken into account [7]. Table 5 
shows results for the process e + e~ — ► 4 jets (ne- 
glecting fragmentation effects) under inclusion of 
several QED and QCD contributions at LEP II 
energies. The input parameters are as in section 
3 and a s = 0.103. The cross section is lowered 
by the initial state QED radiation (ISR), while all 



the other contributions tend to raise it back to its 
Born value. The electroweak (EW) background 
is at the per cent level as well as the interfer- 
ing QCD background (four quark production via 
gluons). In the last row also the non-interfering 
QCD cross section e + e~ — ► 2 gluons + 2 quarks 
is included; it increases the cross section by 16% 
at v^= 175 GeV. 

In figure 1 the distribution of the invariant 



l — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — r 




100 200 

m(i,j) (GeV) 

Figure 1 

Distribution of the invariant masses for e + e~ — 
4 jets at y^s = 175 GeV. 
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masses for the WW signal, the fully corrected 
cross section and the non-interfering QCD con- 
tribution are shown at \fs= 175 GeV. The curves 
have been obtained by fitting histograms from 
EXCALIBUR. The combinatory background is in- 
cluded too, by taking into account any possible 
invariant mass m(ij), so that the area under each 
curve is 6 times the value of the corresponding 
cross section. ISR lowers the distributions and 
ggqq gives a rougly constant positive contribution 
between 30 and 100 GeV. Below the peak the lat- 
ter effect is dominant, while at the peak the sit- 
uation is reversed. A more detailed analysis can 
be found elsewhere [7]. Here, we only report that 
the maximum and the width are not affected by 
QED initial state radiation. On the other hand, 
the QCD and combinatory background increase 
the width. 

6. Conclusions 

Four fermion production in e + e~ annihilation 
will allow important tests of the Standard Model 
and measures at LEP II and beyond. 

There are processes and experimental setups 
for which taking into account only subclasses of 
Feynman diagrams is a good approximation, but 
this is not always the case. As a consequence, 
programs including all diagrams and being able to 
compute any possible four fermion final state are 
indispensable. Also initial state QED radiation 
and QCD contributions turn out to be important, 
and must be taken into account. EXCALIBUR is a 
Monte Carlo program satisfying all these require- 
ments. A Monte Carlo approach is unavoidable, 
since important effects (e.g. the average energy 
loss) are very sensitive to the imposed experimen- 
tal cuts. 
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